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ABSTRACT 

The Ultraviolet Coronagraph Spectrometer (UVCS) observed Doppler shifted 
material of a partial Halo Coronal Mass Ejection (CME) on December 13 2001. 
The observed ratio of [O V]/0 V] is a reliable density diagnostic important for 
assessing the state of the plasma. Earlier UVCS observations of CMEs found 
evidence that the ejected plasma is heated long after the eruption. We have 
investigated the heating rates, which represent a significant fraction of the CME 
energy budget. The parameterized heating and radiative and adiabatic cooling 
have been used to evaluate the temperature evolution of the CME material with a 
time dependent ionization state model. The functional form of a flux rope model 
for interplanetary magnetic clouds was also used to parameterize the heating. 
We find that continuous heating is required to match the UVCS observations. 
To match the O Vl-bright knots, a higher heating rate is required such that 
the heating energy is greater than the kinetic energy. The temperatures for the 
knots bright in Lja and C III emission indicate that smaller heating rates are 
required for those regions. In the context of the flux rope model, about 75% of 
the magnetic energy must go into heat in order to match the O VI observations. 
We derive tighter constraints on the heating than earlier analyses, and we show 
that thermal conduction with the Spitzer conductivity is not sufficient to account 
for the heating at large heights. 

Subject headings: Sun: coronal mass ejections(CMEs) — Sun: activity — Sun: 
corona — Sun: UV radiation 
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Introduction 



The lives on board the Solar and Heliospheric Observatory (SOHO) provides unique 
spectroscopic dia gnostics of CME material through both Doppler shifts and hne intensities 
( iRaymond 1 12002| ) . For CME studies, major questions are the three dimensional structure 
and the energy budget of the ejected CME plasma. UVCS observations contribute to the 
understanding of these two major problems. The temperature evolut ion and energy bal 



2001 



ance during the ejection of CME have been studied by several authors (lAkmal et al. 
Kurnar and Rust1ll996 : Emslie et al. 112004 : Vourlidas et al. [[2000 : Subramanian and Vourlidas 
20071 ). Many studies of CME three dimensional structure have been based on the two di- 
mensional projection of the 3D structure on the plane of the sky and focused on the efforts 
to overcome that loss of information (see Burkepile et al. 2004 and the references in that pa- 
per). In this paper, we present the three dimensional structure, heating, and energy balances 
for a CME observed on Dec. 13, 2001. 

The kinetic properties of CMEs from measurements of white light corona l observations 
have been determined, but suffer from the inaccuracies due to projection effects (IBurkepile et al. 
20041 ). Those make it difficult to study of the three-dimensional structure . To overcorne the 



projection effect , several geometrical assumpti ons such as the cone model (IZhao et al. 



1984 



Zie et al. I l2004l: lYeh. Ding, and Chen 1 120051 ). ice-cream cone model ( iFisher and Munro 



2002 



Xue et al. I l2005l ). and the time difference b etween the appearances of the halo at 



two opposite position angles (IMichalek et al. II2003I ) have been applied to evaluate the CME 



properties. The images reconstructed through polarization analysis indicate that CME is a 
bubble- like structure rather than an expanding lo op arcade structure (ICrifo. Picat and Cailloux 
19831 : iMoran and Davilalbooiloere et al. mO^ . Howard et al. (2007) used SMEI data with 
the assumption of purely radial expansion from the eruption site to estimate the position of 
the ejecta in three dimensions. The three dimensional structure reconstruct ed with UV emis- 



sion lines observed by UVCS shows a ribbon-like structure in one event ( Ciaravella et aL 



20031 ) and a halo CME structure in the view from the solar west in another (ILee et al. II2006I ) . 
The flux-rope model (IKrall and St. Cyr II2006I ) compared statistically with CME ob served 



by So lar Maximum Mission (SMM) shows that the morphology of the CME is hollow (IKrall 



20071 ). The Solar Terrestrial Relations Observatory (STEREO) was launched on 25th Oct. 
2006, and two spacecraft traili ng and leading Earth will provide future studies of the three 
dimensional structure of CME (IKaiser II2005I ). 



The line-of-sight velocity component observed by UVCS makes it possible to estimate 
the angle with the plane of the sky. For evaluating the required heating for the UVCS 
observations, this angle is used to estimate the radial velocities of the ejected plasma, and the 
column density from Large Angle and Spectrometric Coronagraph (LASCO) measurements. 
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The angles with the plane of the sky have been evaluated fo r the 22 halo and parti al halo 
CME cores observed by UVCS to estimate projection effects (ICiaravella et al. II2006I ). 



The energy partitions for flare-CME events show that the CME has the dominant com- 
ponent of the rel eased energy, and t hat it contains a substantial fraction of the available 
magnetic energy (lEmslie et al. II2005I ). Studies of flux rope C ME kinetic energies also show 



that i nternal magnetic energy is a viable source of CME energy (ISubramanian and Vourlidas 



20071 : IVourlidas et al. Il2000l ). Earlier UVCS ob servations of CME found evidence that th e 
ejected plasma is heated long after the eruption (lAkmal et al. 11200 ItlCiaravella et al. Il200ll ). 
and that the total heat going into the ejected plasma is com parable to the ki netic energy. 
Other studies based on t he thermal energy evolution of I CMEs (iLiu et al. II2006I ). emission in 
the EUV bands of EIT (jPilippoy and Koutchmvlbood ) and the ionization states measured 



in ICMEs (IRakowski et al. 



20071 ) also found that CME plasma is strongly heated even after 



it leaves the solar surface. 

In this paper, we find the heating by a procedure similar to that of Akmal et al. (2001). 
We first generate a large grid of models with different initial conditions and forms of the 
heating function. Then, for individual knots observed by UVCS, we select all the models 
that match the observed spectra, and the range of heat inputs for those models is the range 
of heating allowed. 

We find that continuous heating is required to match the UVCS observations. The 
temperature evolution shows a rapid decrease at lower heights and earlier stages, and then 
heating is required to increase the temperature to match the observations. To match the O 
VI bright knots, a higher heating rate is required such that the heating energy is greater than 
the kinetic energy. The temperatures for the O VI, Lya and C III emission indicate that 
different heating rates are required for knots bright in those lines. The required heating rates 
are much larger than wave heating rates for the quiet Sun, and they appear to be larger than 
predicted by thermal conduction. Magnetic heating should occur as the expanding magnetic 
flux rope dissipates energy to reach the simple configuration seen in Interplanetary CMEs. 
The flux rope model of Kumar and Rust (1996) predicts that somewhat more energy goes 
into heat than into kinetic energy. Our results for the knots bright in O VI are compatible 
with that prediction, while the heating rates for the knots bright in Lya are smaller. 

In §2 we describe UVCS, LASCO, Extreme Ultraviolet Imaging Telescope (EIT), and 
Transition Region and Coronal Explorer (TRACE) observations. In §3 we explain the ex- 
pansion model for reconstructing the three dimensional structure from a time sequence of 
long slit spectra and show the three dimensional structure of O VI, Lya, and C III emissions 
observed by UVCS. In §4 we describe the observational constraints and the heating model. 
In §5, the results are presented in terms of the temperature evolution, heating energy, and 
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energy balances. In §6 we discuss our results and compare the heating rates with those 
expected for thermal conduction. In §7 we summarize our studies. 



Observations 



UVCS ( iKohl et al. Ill995l ) obtains spectra of the solar corona inside an instantaneous 
field of view (FOV) given by the 40' long spectrometer entrance slits and the slit width (21" 
in this observation), which can be placed between 1.5 and 10 -Rq. The UVCS slit was pointed 
above NOAA 9973 (N16 E09 at 14:20 UT on Dec. 13 2001) at a position angle P.A.=349° at 
2.4 Rq. The P.A. and height are the values at the point along the slit that is closest to Sun 
center, which is 5/9 of the way from the east edge of the slit. The series of 300s exposures 
separated by about 30s readout time began many hours before the event and lasted until 
19:16 UT at the same position. The data were obtained with a spatial binning of 6 pixels 
(42") and a spectral binning of 1 pixel, 0.0993 A and 0.0915 A for primary and redundant 
wavelengths, respectively. 

Figure 1 shows the composite image of EIT 195, UVCS, and LASCO C2 observations. 
A bright flare was observed in the EIT 195 A band beginning at 14:24 UT, and the peak 
intensity was seen at 14:36 UT in the 12 minute cadence images. In this event, TRACE 
observed this active region in the UV continuum (1600 A and 1700 A) and C IV 1550 A. A 
filament eruption was observed in the UV continuum, 1600 A around 14:20 UT (Figure 2). 
The leading edge of the CME first appeared in the LASCO C2 FOV at 14:54 UT, and in 
the C3 FOV at 16:18 UT. The Geostationary Operational Environmental Satellite (GOES) 
X6.2 class flare began at 14:20 UT, peaked at 14:30 UT, and lasted until 14:35 UT. Figure 
3 shows the observation time of each instrument. EIT observed this active region with 12 
minute cadence. LASCO C2 and C3 obtain images with 24 minute cadence, however the C2 
and C3 data were not available from 14:54 UT - 16:30 UT and 14:42 UT - 16:18 UT during 
the CME expansion, respectively. At 14:54 UT, C2 observed the leading edge of this event 
and, then at the next observation, the CME has passed beyond of C2 FOV. 

From 14:56 UT, UVCS observed Doppler shifted material in the O VI 1032 A, Lya 1216 
A and C III 977 A lines. Also the faint spectral lines, O V] 1218.35 A, and [O V] 1213.9 
A, N III 989.90 A, N III 991.58 A, C II 1036.34 A, and C II 1037.02 A were observed in 
several exposures. It is especially important that the observed O V] (2s^ ^Sq — 2s2p ^Pi) 
intercombination line (violating the selection rule A S=0, notation: right square bracket, O 
V]) at 1218.35 A and the [O V] (2s^ ^5*0 — 252^ ^P2) forbidden line (notation: square brackets, 
[O V]) at 1213.9 A can be used to determine the plasma density (see §4.2.1). Because the 
Einstein A value of the forbidden line is 0.022 s~^, the line is quenched at densities above 
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about 10^ cm~^ as shown in Figure 9 of Akmal et al. (2001). To investigate the heating 
rates, we analyzed the spectral lines for the fourteen blobs observed in four exposures (refer 
to Table 1). Figure 4 shows the UVCS observations of O VI (each left panel) and Lya, O 
V], [O V], N III (each middle panel), and C III (each right panel) from 14:56 UT to 15:12 
UT. The Lya and O V lines are observed in the redundant wavelength with the N III line 
in the primary wavelength. The positions along the slit of the fourteen blobs are marked at 
the right axis of the middle panel with white bars. The intensities and line-of-sight velocites 
(Vlos) of each blobs are presented in Table 2 (O VI and O V) and Table 3 (Ly«, C III, 
C II, N III). For instance, in the case of blob A at the first detection of the CME at 14:56 
UT, we selected the Doppler shifted O VI emission with -0.98 A and -1.68 A (evaluated by 
two Gaussian fit) indicating line-of-sight velocities of 286 km/s and 459 km/s, respectively 
at P.A. = 361 ° - 366°. In the exposure at 15:01 UT, we selected three blobs (B.l, B.2, and 
B.3) and those are also marked on the left panel of Figure 4. The C III emission of blobs A 
and B was detected at the edge of the observed wavelength range in this data panel. In the 
case of region C observed at 15:07 UT, we selected six blobs, one detected only in O VI and 
C III (C.l), one detected in the bright emission of Lya but relatively faint emission of O VI 
(C.2), and others detected in the emission of the forbidden [O V] A1213.9 line (C.3 - C.6). In 
the 15:12 UT observations, we selected four blobs (D.l - D.4) detected also in the emission 
of [O V] . The line intensities for most of blobs were evaluated by fitting two Gaussians to 
the line profile. 



3. Three-dimensional reconstruction 

A three-dimensional reconstruction of the 13 December 2001 CME has been made based 
on O VI 1032 A, Lya 1216 A, and C III 977 A observed from 14:56 UT to 19:16 UT. Briefiy, 
for each pixel along the UVCS slit, we first determine the velocity components in the plane 
of the sky for each exposure by dividing the difference between the eruption time and the 
time of the exposure into the difference in position between the pixel and the location of 
the eruption. The observed Doppler shift provides the line of sight velocity component 
for each pixel in each exposure. When different structures are superposed along the line 
of sight, they may show up as separate velocity components and can thus be separated 



into, for instance, structures on the front and back sides of the expanding CME (ILee et al. 



20061 ) . which were separated with multiple Gaussian fits. We then use the 2 plane-of-the- 
sky velocity components and the line-of-sight velocity to project the structure forward in 
time, adding in the material observed in each successive exposure. For this event, due to 
the 300 s exposure time, blobs with different properties were sometimes superposed in the 
same pixel during the same exposure. The velocities in the plane of the sky would cause 
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those blobs to be superposed at later times in the 2D projection onto the plane of the sky, 
but the different line-of-sight components of the blobs make it possible to separate them in 
the 3D structure. Our method essentially assumes that the expansion is self-similar. This 
assumption is checked by comparing the projection of the 3D structure onto the plane of 
the sky to LASCO images at later times. We use the same assumptions as the recent 3D 
reconstruction of UV lines for the 2002, April 21 CME observed by UVCS by Lee et al. 
(2006) who give a detailed description of the method. This model has two key assumptions. 
One is homologous expansion in all directions. The other is that the line-of-sight velocities 
are only caused by the material's expansion, rather than thermal or turbulent motions. 



3.1. Pre-processing of the data 



We subtracted as background the average of 35 pre-CME exposures observed from 
08:02 UT to 14:16 UT. The scattered Lya continuum was eliminated by a linear equation 
for the background emission in the observed wavelength ranges. T he Gaussian wings on the 
instrumental line profile were corrected by an iterative procedure ( iKohl et al. 1119971 ). Since 
the detector has some distortion, we evaluated the reference line center of the O VI and 
Lya by Gaussian fitting after binning by three spatial bins along the slit using the pre-CME 
observations above. The line center of the C III was evaluated using the O VI line center. 



3.2. Expansion Model 

The onset position of the CME eruption is obtained from the position of NOAA 9973, 
which was located at 0.084 Rq to the west (hereafter the x-coordinate) and 0.218 Rq to the 
north (hereafter the z-coordinate) according to the EIT image at 14:24 UT. The position in 
the line-of-sight direction was thus -0.972 Rq toward the Earth (hereafter the y-coordinate) 
assuming the Sun is a sphere {Rq=1). First, We assume that each element of plasma moved 
in a straight line from the eruption site to the UVCS pixel where it was observed, as described 
above. Second, we assume that all material erupts at the same time and location, and it 
has constant speed before reaching the UVCS slit at each spatial bin. This implies that the 
plasma reaching the slit during later exposures travels more slowly in the plane of the sky 
than the plasma that reaches the slit earlier. The ejection is spread out in both location 
and time, but only by about 100 arcseconds and 10 minutes in time as estimated from the 
TRACE movie. Given the limitations of spatial resolution and cadence of the UVCS and 
LASCO data, the assumption of a single point in space and time is adequate. 
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3.3. Images from the three-dimensional reconstruction 

A movie has been made using a 40 Rq cube with 0.1 Rq pixels (see electronic version). 
The reconstructed movie shows the distributions of the O VI, Lya, and C 111 emissions 
(Figure 5). The red, green, and blue colors correspond to the O VI, Lya, and C 111 emission, 
respectively. The yellow (O VI and Lya), sky (Lya and C 111), violet (O VI and C 111), 
and white colors (O VI, Lya, and C 111) show the overapping features of those lines. At the 
beginning of the observation, 14:56 UT and 15:01 UT, the blue shifted components of C 111 
emission were detected at the end of the wavelength range of this data panel (marked in 
Table 3). This could underestimate the line-of-sight velocity to the Earth in the reconstructed 
movie. 

We present four views from the solar east-Earth (upper left), the solar east-backside 
of the Sun (upper right), the backside of the Sun (lower left), and the solar west (lower 
right). The reconstructed structure is extended towards the Earth, as indicated by the 
strong blue shifts. The views from the solar west (lower right) make it easy to understand 
the distribution of the three emission lines along the line of sight. In the CME front, the O 
VI and Lya emissions are positioned close together (red, green, and a mixed yellow colors). 
Following those, the Lya and C 111 emissions are closer to the middle of the ejected CME 
plasma (green, blue, and a mixed sky colors). The C 111 emission is seen in the innermost 
part of t he reconstructed mo vie, in agreement with earlier indications of a core of cooler 



material (lAkmal et al. Il200ll ). The view from the backside of the Sun to the Earth (lower 
left panel in Figure 5) also shows the C 111 emission in the innermost part of the ejected 
CME material. The bright features of O VI, Lya, and Clll observed by UVCS show the arch 
(shell-like) structure in the view from the solar backside (please see the movie and the image 
at the bottom left in the Fig. 5). Inside of the shell-like structure, it is much fainter and 
seems to be void. Kr all (2007) compared the C MEs observed by SMM with a parameterized 



3D flux rope model (IKrall and St. Cyr II2006I ) and found that the CMEs observed by SMM 



in his study are consistent with the hollow flux-rope geometry. Our 3D reconstruction looks 
like the hollow flux rope structure, as best seen in the view from the back of the Sun. 



4. Energy Budget 

To investigate the energy balance of the 13 December 2001 CME plasma, we used the 
fourteen blobs from the 4 exposures observed by UVCS from 14:56:19 UT to 15:12:48 UT 
(Table 1, see also Figure 4). The bright features of the fourteen blobs occupy small spatial 
regions along the slit in one or more of the spectral lines. We analyze all the observed lines 
along the bright features which are marked with the white bar at the edge of the middle 
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panel in Figure 4. 

The procedure is much hke that of Akmal et al. (2001), but for this event we have 
different constraints because different hnes were observed, and because we use the three- 
dimensional reconstruction to constrain the depth along the line of sight. First, we evaluate 
the expansion of the CME plasma based on the TRACE, LASCO, and UVCS observations 
(§4.1). Second, we compute grids of models for a broad range of initial temperatures, den- 
sities, and forms of the heating function to calculate the ionization fraction along the CME 
expansion. The ionization fraction is evaluated by a time-dependent ionization code using 
density and temperature as input which arc calculated from the expansion law (§4.2.2) and 
the energy equation (Eq. (5) in §4.3). The ionization fractions are also used to compute 
the radiative cooling rate. At the final step of each model, the density, temperature, and- 
ionization fraction are used to evaluate the line intensities (Eq. (2)). Finally, the line inten- 
sities evaluated from each heating model are compared with the O VI, O V], Lya, C III, C 
II, and N III line emission observed by UVCS. 

The heating laws for fast solar wind (Allen et al. 1998), in which proton and electron 
rates differ, as well as a density-proportional heating rate are used to model the heating rates. 
The flux rope model for interplanetary magnetic clouds (Kumar and Rust 1996) is also used 
to evaluate the energy balances of the ejected plasma. The heating rates in the models 
contain all forms of thermal energy input including wave dissipation, magnetic heating and 
the divergence of thermal conduction, which we will consider in Section 6. Here we model 
the ejected plasma including radiative and adiabatic cooling to investigate the temperature 
evolution of the CME material with a time dependent ionization state model. We describe 
the model in more detail in the following subsections. 

4.1. Expansion Model 

We use an acceleration profile to evaluate the height of the ejected plasma rather than 
the constant velocity that was used to make the three-dimensional reconstruction. Chen et 
al. (2006) showed that the distance of the two footpoints in the active region is related to 
the maximum expansion velocity of the CME flux rope. We determined the peak time of the 
acceleration, 14:25:00 UT, based on the distance of the two footpoints from TRACE 1600 
A observation. Then the start time of the expansion, 14:10:12 UT, was evaluated to match 
the heights of the CME plasma observed by TRACE, UVCS, and LASCO observations. 
Different slopes of the acceleration were used for each exposure. The applied accelerations 
are 900 m/s^, 765 m/s^ for the A and B blobs, respectively and 630 m/s^ for the C and D 
blobs. The TRACE 1600 A observations also show that the brightening in the active region 
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starts around 14:10 UT and increases slowly after that. The evaluated heights of the ejected 
CME plasma are presented in Figure 6. 



4.2. Observational Constraints 

Electron density, kinetic temperature, and line intensity are used to constrain the ac- 
ceptable models for a broad range of initial temperature, density, and three heating functions 
for each of the fourteen blobs. The electron density evaluated at the UVCS slit height is 
used to determine the electron density at each time in the expansion model and to evaluate 
the emission measure (EM). The kinetic temperature of Lya is used as the upper limit to 
the temperature at the UVCS slit height. The observed line intensities are compared with 
the intensities obtained from the models. 



4-2.1. Electron density 



The ratio of the [O V]/0 V] lines is a reliable density diagnostic. This ratio has been 
used to determine the electr on density of the ej ected CME plasma, and it is most useful 
for densities of 10^-10^ cm'^ flAkmal et al. IboOlh . In this event, UVCS observed the [O V] 
forbidden line at 1213.85 A and the O V] intercombination line at 1218.39 A in several 
blobs, though the 1213.85 A line could not be cleanly separated from the N III 991.85 A 
line in most cases. For the seven blobs observed at 15:07 UT and 15:12 UT, the electron 
number densities were evaluated from [O V]/0 V]. We present the line profiles of the three 
blobs C.5, C.6, and D.3 in Figure 7. Both primary and redundant wavelengths are presented 
together. The [O V] line may be contaminated by N III 991.58 A, and while we attempted to 
remove the N III by assuming that its profile is the same as that of C III 977 A, we consider 
the [O V] intensities to be upper limits and the density estimates to be lower limits. We 
use the electron number densities from the ratio of the [O V]/0 V] for the three blobs, C.5, 
C.6, and D.3, for the anlysis of the heating rates and those are marked with '^^ in Table 2. 
The electron density eval uated from the O V ratio and version 5.2 of CHIANTI (jPere et al. 
I997I : Ihandi et al.lbooeh is presented with uncertainties in Figure 8. For other blobs, the 
electron number densities are evaluated from the column densities observed by LASCO, 
divided by the line-of-sight depth based on the O VI line emission. The three dimensional 
reconstruction of the O VI lines shows the reconstructed width is similar to the CME width 



from the LASCO observation (ILee et al. II2OO6I ). We take the depth of the O VI line to be 
the FWHM multiplied by the travel time from the source region to the UVCS slit. The 
evaluated electron number density and column density were presented in Table 2. 
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The LASCO C3 observation at 16:18 UT was used to evaluate the column density along 
the line-of- sight. Unfortunately LASCO didn't observe this event when the [O V] and O V] 
lines were observed from 15:07 UT to 15:12:48 UT (refer to Figure 3). The three-dimensional 
reconstruction makes it possible to estimate the height and the angle with the plane of the 
sky. The heights of the UVCS blobs predicted by the constant velocity reconstruction are 
higher than their positions in the LASCO C3 image at 16:18. We use the positions from the 
reconstruction using the observed acceleration to connect the UVCS blobs to the LASCO 
image. In Figure 9, the reconstructed positions of the observed blobs are presented with 
the LASCO C3 observation, and each box represents the A, B, C, and D blobs in the Table 
1. The reconstructed positions of blobs, A— D are also indicated by '3D position' in Figure 
6. Vourlidas et al. (2000) showed how the densities derived from LASCO increase with 
the angle with the plane of the sky due to the angular dependence of Thompson scattering. 
Plane-of-thc-sky angles of 40° ~ 50° were found for the observed blobs from the line-of-sight 
velocities. We used the angle with the plane of the sky of 45 ° to evaluate the column 
densities from the LASCO observation. The column densities were calculated by averaging 
the area of each box in Fig 9. After the column density was obtained, we scaled the density 
inversely with height squared to estimate the column density at the height of the UVCS 
observations. 



4-2.2. Temperature and line intensity 

The line widths obtained by the Gaussian fit for the fourteen blobs contain the kinetic 
temperature and non-thermal (bulk-velocity) components. The upper limits on the kinetic 
temperatures of O VI and Lya are presented in Table 2 and Table 3, respectively. Note that 
in the 3-D reconstruction, the line width is assumed to be caused by the material's expansion 
along the line of sight. Here we regard the kinetic temperature to be the upper limit of the 
electron temperature at the UVCS slit height as one of the criteria for selecting the acceptable 
models. The Lya width constrains only the ion temperature in the Lya emitting region since 
the O VI and Lya generally do come from different regions. The temperature Tj- is defined 
as: 

n = Ijvy^' (1) 

where m is the mass of the ion and vi/e is the velocity derived from Doppler half width, 
AAi/e. We subtract about 0.3 A (3 pixels) in quadrature to account for the effect of the slit 
width of 0.298 A (0.275A for the redundant channel). 
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Radiative and adiabatic cooling along with the heating rates are used to investigate the 
temperature evolution of the CME material in the time dependent ionization state model. 
We use the expansion law, riuvcs/inQ — riuvcs) oc (t/tuvcs)"', which determines the adiabatic 
cooling rate {Ladiabatic in Eq- (5)) from T /n^~^ (Akmal et al. 2001, refer §4 in their paper). 
Initial densities (no) are examined from 3 x lO^cm"^ to 2 x lO^^cm"^ with initial temperatures 
from 1 X 10^ to 6.3 x 10^. tuvcs is the duration time of the CME expansion from 14:10:12 UT 
(refer §4.1) to the UVCS observations. Using the density ijiuvcs) determined at the UVCS slit, 
2.4 i?0, the expansion is assumed to be a power law. We adopt the power law of a=3, which 
means an expansion in both length and radius and a slower expansion at the earlier stages 
than the expansion with a=l or 2. Since we don't know exactly how the CME has expanded 
from the origin to the solar corona, the power law assumption which has been adopted in 
the earlier work (Akmal et al. 2001) is the simplest way to parameterize the density falloff 
with height. Also the choice of the a=3 matches the observed slower expansion in the earlier 
stages. During the expansion of the CME plasma, the ionization fraction is evaluated from 
the time-dependent ionization state model. The ionization and recombination rates were 
taken from an updated version of the code of Raymond (1979). This atomic physics package 
computes the coUisional ionization, radiative recombination and dielectronic recombination 
rates for the astrophysically abundant elements, and the evolution of the ionization states of 
the elements was computed with the from the usual set of equations for the rate of change of 
each ionic fra ction. The time-dependent ionization balance was used with solar photospheric 



abundances (iGrevesse and Sauval Ill998l ) to compute the radiative cooling rate (L^, in Eq. 
(5)). Photospheric abundances were chosen because the knots we analyze are probably 
ejected prominence material, rather than coronal plasma. The computed ionization fraction 
is used to determine the line intensity. In coronal conditions with the isothermal assumption, 
the line intensity is expressed as follows (e.g. Ko et al. 2006). 



Hine 



— —I e{T^)dEM{Te) photons cm'^ s'^ sr"^ (2) 
An nn J 



Where n^i/nH is the abundance of the element relative to hydrogen. In our model, the line 
intensity is calculated for a single temperature at any time. So we use the emission measure, 
EM (= 0.8?T,eA'"e) instead of dEM without the integral. 0.8 is the hydrogen abundance 
relative to the electron density, rtf. is the electron number density determined from the O 
V lines or from LASCO observation divided by line-of-sight depth (see §4.2.1), and N^. is 
the column density determined from the LASCO image (see §4.2.2). e(Te) is the emissivity, 
which is defined as 

n ' 

e(Te) = -^{Te)qii„e{Te), (3) 
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nion/nei is the ionization fraction, qiine is the colhsional electron excitation rate given 

by 

y J- e 

where fi^ is the coUision strength taken from version 5.2 of CHIANTI and Wi is the 
statistical weight. 

The observed line intensities are compared with the line intensities predicted from the 
three heating models (see §4.3). To match the O VI blobs, we use the ratio of the O VI/ 
O V] which allows a factor of two margin because the density obtained from O V may 
not pertain to the O VI emitting gas. To match the Lya, C III, C II, and N IIII blobs, we 
require that the predicted emission is at least as large as the observed values, but it is smaller 
than three times the observed intensity. We allow this margin because of the difficulty in 
evaluating the emission measure, given that we do not have a density from the O V lines 
for this cooler gas, and the combination of line-of-sight depth with LASCO column density 
carries a greater uncertainty. When the C III was observed near the edge of the panel, we 
exclude this constraint. In cases where the O VI/0 V] ratio is not well matched (C.4 and 
D.4 for the heating model of the Kumar and Rust, marked by in Table 4), we use the same 
constraints for the Lya, C III, C II, and N III blobs. We present the detailed constraints 
and the evaluated energy in Table 4. 



4.3. Heating model 

Earlier UVCS obser vations of CMEs found evidence that the eje cted plasma is heated 



long after the eruption (lAkmal et al. I l200ll : ICiaravella et al. I l200ll ) . We investigate the 



heating rate for the following three heating models (refer to §4.3.1 and §4.3.2). The heating 
rate, H, enters the energy equation. 



5 dT _ 

2^ at ~ ~^'^^P^^ ~'~ ^a.d,iabatic ~ ^) 



(5) 



where n is the density of both protons and electrons, k is the Boltzmann constant and and 
Up are an electron and proton density, respectivly. L,. and Ladiabatic are radiative cooling rate 
and adiabatic cooling rate, respectively. H includes all contributions to the heating, such as 
wave dissipation, magnetic reconnection, turbulent heatin g or the divergence of conductive 
flux, or shock waves generated by the reconnection outflow (jShiota et al. II2005I ). We combine 
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turbulent heating with wave dissipation, so 



wave 



(6) 



We have used the physical value of the specific heat ratio 7=5/3. Smaller values of 7 
are often used in coronal models with the heating rate set equal to zero in order to obtain 
plausible temperature profiles without knowledge of the heating rate. In this case we aim to 
derive the heating rate, so the value of 7 for a monatomic gas is appropriate. 

We use a wide range of heating rates with three forms of the heating function. Two of 
them are taken from parameterized heating rates for the solar corona, and the third is from 
the CME model of Kumar and Rust. 



The heating law for fast solar wind (Allen, Habbal and Hu, 1998), in which proton and 
electron heating rates differ (hereafter denoted as Q oc Qahh), and the density- proportional 
heating (hereafter Q oc n) are used to evaluate the heating rates. 

The density proportional model was chosen, as in Akmal et al. (2001), to provide a 
heating rate that declines more gradually than the exponential rate of Allen et al. (1998). 
It implies no specific physical mechanism, but for instance saturated thermal conduction 
with a constant temperature gradient might give such heating. Shiota et al. (2005) have 
described how the impact of reconnection outflows on the trailing edge of a CME core can 
generate slow mode shock waves that travel around and into the core. Lacking knowledge of 
the evolution of the reconnection outflow or the location of the shock dissipation, we cannot 
predict the functional form of this heating rate, but it is plausible that it declines on a time 
scale comparable to the travel time to the height of the UVCS slit. 

For the model, Q oc Qahh, the proton and electron heating rates (Qp and Qg) are 
defined as. 



Qe = Qeoe-('-^^)/"^ ,Qp = Qpoe-^'^-'''^/''^+qOO{ — ] (1 - e-^^-^^V^^^), (7) 



in their paper. We consider a constant heating rate with an exponential falloff with scale 
height 0.7 Rq ((7e,p). Qeo, Qpo, and qOO are the strength of heating terms taken from the 
SW3 model in their paper, r is the distance from the initiation site and R5 is solar radius. 



4.3.1. Parameterized heating 
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The temperature changes due to the proton and electron interaction are considered for the 
heating model, Q oc Qahh- We model the heating functions for a range of constants of 
proportionality covering three orders of magnitude. 



4-3.2. Magnetic Heating 

The magnetic clouds in Interplanetary Coronal Mass Ejections generally have simple 
magnetic structures that can be fit with a Lundquist solution (e.g. Lynch et al. 2005), 
a minimum energy state for a given magnetic helicity. To reach this state, the complex, 
stressed magnetic configuration of the erupting structure must relax by dissipating magnetic 
energy. The MHD model of Lynch et al. (2004) shows that this occurs by the time the CME 
reaches about 15 Rq. MHD models do not generally show predicted heating rates, and in 
any case the heating is unlikely to be uniform throughout the CME structure. However, the 
model of Lynch et al. (2008) indicates that about 15% of the total change in magnetic free 
energy goes into kinetic energy. 

In order to parameterize magnetic heating we use the flux rope model for interplanetary 
magnetic clouds of Kumar and Rust (1996; hereafter Q cx Qkr)- The model predicts 
that the magnetic energy goes into overcoming solar gravity, kinetic energy for the CME 
expansion and heating energy. An important feature of the model is its use of the principle 
of conservation of magnetic helicity. It suggests that the magnetic energy stored in an 
expanding plasma should decrease with expansion, {Hm oc length scale x Um), where Hm 
is the magnetic helicity and Um is the magnetic energy. Kumar and Rust (1996) assume 
that the flux rope evolves through a self-similar series of axisymmetric states each having 
the lowest energy consistent with conservation of magnetic helicity. The assumption of 
self-similar axisym metric expansion rn ay be adequate when the CME is much larger than 



the source region (jUralov et al. 1120051 ). but it is not likely to be accurate at small heights 



(lAkmal et al. Il200ll ). Kumar and Rust (1996) also assume that magnetic energy is dissipated 



in the expansion process, though they do not identify the dissipation mechanism in detail. 

We use the length scale, / = 27ra, where a is the radius of the flux rope (refer to Figure 
2 in their paper). We use the distance of the CME plasma from the initiation site as the 
radius 'a'. Kumar and Rust found that the change of magnetic energy is dUm = Umilollodl/t^ 
and that as the CME expands the magnetic energy lost is partitioned between kinetic + 
gravitational and thermal energies according to a parameter s = 1 — sin{6o)/'n' that may 
vary from one event to another. 6o is an effective toroidal angle (refer Fig 5 in their paper). 
As stated in their paper, the kinetic energy of major radial expansion about the center of 
mass is same as the kinetic energy of center of mass (the kinetic energy of the minor radial 
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expansion is negligible for a large aspect ratio toroidal cloud) and the force acting on the 
center of mass is maximum on the leading edge of the flux rope. With twice the center of 
mass kinetic energy derived by integrating of the motion of the center of mass from the initial 
distance to the leading edge of the flux rope {9q ^ vr/2), they presented energy conservation 
as dQ + dUm + dUc + dUKE = 0. Using the energy conservation, we derive the heating energy 
during the expansion of the CME plasma. If dUc is the change in gravitational energy and 
dUxE is the change in kinetic energy, then 

dQ = -hdUm (8) 

where h — (l + dUG/dUm + dUKE/dUm) is the fraction of the lost magnetic energy appearing 
as heat. We find the kinetic and gravitational energy from the velocities and distances 
explained in §4.1. The final kinetic energy and the initial gravitational energy are fixed by 
the observed velocity and an initial height of CME plasma, respectively. We find the heating 
rates to match the observations for in the range of 10—90 °. The blobs we observe appear 
to be part of the ejected prominence, and as such they are part of flux rope. The energy 
budgets we obtain pertain specifically to these blobs, but lacking any reason to assume 
otherwise, we believe that the heating is uniform throughout the flux rope. 

5. Energy balances 

5.1. Temperature evolution 

Figure 10 shows the temperature evolution of both protons and electrons when one of 
the acceptable heating rates is applied for matching the O VI and O V] line intensities of blob 
C.5. The applied heating rate is 63 times that for the fast solar wind model (QocQahh)- The 
upper two panels show the proton and electron temperature evolutions for all applied initial 
densities (3.16x10^—1.99x10^^ cm^^) for the particular initial temperature, 1.58x10^ K. 
The higher initial densities (red and yellow ranges in the rainbow colors) drop more rapidly 
than the lower initial densities, but all initial densities reach the same temperature at the 
time of the UVCS observation. The lower two panels show the same feature as that above 
for all applied initial temperatures (1.x 10^ — 6.3x10^ K) at that particular initial density, 
5x10^ cm~^. The proton temperature increases more than the electron temperature because 
of the different heating rates for the protons and electrons (see §4.3.1). The temperatures 
drop rapidly at the lower heights and in the earlier stages of the expansion, and then they 
remain fairly constant depending on the heating rates. This is a similar result to that of 
Akmal et al. (2001) (see §5.2 and Figure 12 in this paper, also see Figure 12 in their paper). 
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Our models show that the temperatures at coronal heights are determined by the heating 
rates rather than the initial temperatures and densities. In Figure 11, we present the line 
intensities (upper panels) and the line ratios (lower panels) for all initial densities with the 
temperature predicted from a range of heating rates which are 0.63 - 630 times Qe and Qp 
of the QccQahh (initial temperature = 1.58x10^ K). The temperatures that agree with the 
O VI, Lya, C III, C II, and N III emissions indicate that different heating rates are required 
for those ions. For this model, the heating rates for the O VI blobs and for the blobs of 
cooler ions are 63 times and 2.5 times for the fast solar wind model, respectively. Hereafter 
HI represents the heating rates that match the O VI emission and H2 represents the heating 
rates for matching the Lya, C HI, C II, and N HI emission. 



5.2. Heating, cooling, and energy budget 

Figure 12 shows the evolution of the temperature, density, and ionization state for one 
of the acceptable models explained above. The temperature and the ionization fractions 
change rapidly at the lower heights. For the HI models, the O VI and O V ionic fractions 
maintain higher values after the rapid increase at lower heights, but the fractions of other 
ions, Lya, C III, C II, and N III emissions, decrease again (the middle panel in Figure 12). 
For the H2 model, the fractions of the O VI and O V decrease and the others maintain their 
fractions (the right panel in Figure 12). The rapid change in the ionic fration at low heights 
reflects the change in the eletron temperature when they are in ionization equilibrium. On 
the contrary, at heights higher than 2 R©, the rapid expansion tends to 'flatten' the ionic 
fraction indicating that the ions are frozen-in. The outflowing ions flow too fast and the 
density is too low for them to exchange electrons between neighboring ionization states from 
collisions. Figure 13 shows the heating and cooling rates for the HI model in the upper panels. 
The adiabatic cooling rate increases rapidly at the earlier stages (the upper right in Figure 
13). This drastic cooling rate requires that continuous heating increases the temperature to 
match the observations. The thermal and radiative loss energies arc smaller than the heating 
and adiabatic cooling energies. The heating energies for the acceptable models for each of 
the three heating functions are shown in Figure 14. The heating rates for the O VI (HI) and 
for the Lya, C HI, C II, and N HI (H2) observations are presented as diamonds and crosses, 
respectively. HI is over ten times higher than H2 and the QocQkr model does not match 
the observations for the Lya, C HI, C II, and N HI blob. The integrated heating energy 
balances cooling energy, and the total heating is greater than the kinetic energy in both the 
QccQahh and QccQkr models, as shown in the lower right in Figure 13 and Figure 15. The 
Kumar and Rust model predicts that the magnetic energy decreases conserving magnetic 
helicity. 



-17- 



6. Discussion 

In Table 4, we present the heating and thermal energy for all fourteen blobs. The kinetic 
energy of this CME is 1.6-3.3x10^^ erg/g which is evaluated from the velocity of the CME 
plasma at the UVCS slit, 2.4 Rq. The gravitational and ionization energy are 5.6x10^^ 
erg/g and 1.9—2.1x10^^ erg/g, respectively. In cases where the model does not satisfy the 
constraints from the observations, these are marked as a - e (details in Table 4). The 112 
model could not satisfy all the constraints together for the several blobs because the lines 
Lya, C III, C II, N III might come from different regions along the line of sight. For all the 
acceptable models, the heating energy for HI is higher by nearly one order of magnitude 
than that for H2, and it exceeds the kinetic energy. A recent study of heating rates for 
Interplanetary coronal mass ejection (ICME) from Helios 1 and Helios 2 observations found 
that plasma turbulence heats the proton and alpha particles and determined t he heating 
rates to be 1.9 xlO^erg g"^ s~^ at 0.3 AU and 1.9 xlO^erg g-^ s"^ at 20 AU Jliu et al 



20061 ). The heating rate, 1x10^^ erg g^^ s^^ at 0.01 AU, estimated by extrapolating from 
their study (Figure 6, in their paper) is similar to our result of the heating rate, 5x10^^ erg 
cm~^ sec~^ at 2.4 R© (taking the electron density 1.x 10^ cm~^ and the mass from assuming 
10 % Helium, 1.974x10-^4 g). 

The heating energies from the Kumar and Rust model are also greater than the kinetic 
energy. For all acceptable models, 70-90% of the magnetic energy goes into the heating 
energy, in agreement with their theoretical calculation. The December 13 2001 event is 
probably a prominence eruption and it is presumably defined by magnetic field (refer to 
Figure 2). The fractions of the magnetic energy that go into heat for the earlier observed 
blobs are higher than these later ones. For H2, we couldn't find a heating rate that satisfies 
the intensities of those low-temperature lines. It is possible that the CME plasma inside the 
flux rope does not follow the heating law of the Kumar and Rust magnetic cloud model or 
that the assumed self-similar expansion is not a good approximation in the early stages of 
the eruption. 

Figure 16 shows the time scale of radiative cooling, adiabatic cooling, and thermal con- 
duction for one of the acceptable models for HI of C.5 (Q oc n). We can use the temperature 
gradients from this model with the Spitzer thermal conductivity to evaluate the importance 
of the thermal conduction contribution to the heating. Larger temperature gradients would 
give too low a temperature at the height of the UVCS slit, so this gives an upper limit to the 
thermal conduction contribution. Other parameterizations of the heating rate could differ 
slightly, but not at the level needed for thermal conduction to balance adiabatic cooling. 
The thermal conductive flux and time scale can be estimated as 
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/{k+i-U-i) (9) 



and 



respectively. Here i is time step and V-Fc is the divergence of the conductive flux, is 
1.8xlO~^/lnA erg cm~^ K~^/^ sec~^, and In A is Coulomb logarithm. In the earliest stage of 
the expansion, the time scale of thermal conduction is much smaller than other cooling time 
scales. In this stage, the thermal conduction would maintain the initial high temperature 
and might change the temperature structure during the first 200 seconds and flatten out the 
initial steep drop (or rise) in temperature if the temperature gradient lies along the magnetic 
field. In the plot of the thermal conduction time scale, the peak around 300 seconds is due 
to the change of the sign of the divergence of the conductive flux. Those plots indicate that 
the adiabatic cooling is much more important than other terms after earlier stages. We 
have assumed classical Spitzer conductivity in equation (9). Saturated conductivity would 
provide a smaller heating rate, but it is possible that a non-Maxwellian velocity distribution 
gives a larger conductive heating rate. 

Due to the expansion of CME, the plasma cools rapidly at earlier stages, and then 
the heating energy becomes the major contribution in the energy budget of the ejected 
CME plasma, even greater than the kinetic energy. Guhathakurta et al. (2006) studied 
the heating of the corona and they also found that the thermal and non-thermal heating 
terms contribute to the temperature and heat flux in the low corona (refer Figure 4 in their 
paper). A similar result for the heating energy of CME plasma based on the Advanced 
Composition Explorer (AC E) observations shows that the plasma requires further heating 



following filament eruption (IRakowski et al. 1120071 ) 



Our results confirm the studies of Akmal et al. (2001) for a different event and obtain 
a narrower range of the heating rates in both ions and electrons. Also our study shows that 
heating rates are consistent with Kumar and Rust's model. About 75 % of the magnetic 
energy goes into the heating energy for the O VI observation with Kumar and Rust's model, 
in agreement with their prediction. 
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7. Summciry 

UVCS observed Doppler shifted material in the O VI 1032 A, Lja 1216 A, and C III 
977 A hnes on December 13 2001. Fainter spectral lines, O V] 1218.35 A, [O V] 1213.9 A, 
N III 991.58 A, C II 1037.02 A, were observed in several exposures. A three-dimensional 
reconstruction of the 2001 December 13 partial halo CME has been made based on the O 
VI, Lya, and C III lines. The reconstructed structure shows the distribution of those ions 
as viewed from all directions. The structure viewed from the backside of the Sun towards 
the Earth looks like a hollow flux rope. The C III emission is located in the innermost part 
of the ejected CME material. The reconstructed structure is also used to find the column 
density with LASCO observations. We investigate the heating rates by a procedure similar 
to that of Akmal et al. (2001). For the individual knots observed by UVCS, we generate a 
large grid of models having different initial conditions and forms of the heating function. We 
find that continuous heating is required to match the UVCS observations. The temperature 
evolution shows a rapid decrease at lower heights and earlier stages, and then heating is 
required to increase the temperature to match the observations. To match the O VI bright 
knots, a higher heating rate is required such that the heating energy is greater than the 
kinetic energy. The temperatures for the O VI and Lya, C III. C II, N III emission indicate 
that different heating rates are required for the bright knots in those lines. About 75 % of 
the magnetic energy goes into the heating energy for the O VI observation with Kumar and 
Rust's model, in agreement with their prediction, but those models do not match the cool 
line blobs. 
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Table 1. Blobs observed by UVCS 



Blob No. 


UT 


P.A. 






Observed Lines 


A 


14;56;19 


361°-366° 


o 


VI, 


Lya, C III, C II A1036.34, C II A1037.02, N III A991.58 


B.l 


15:01:47 


363°-364° 


O 


VI, 


Lya, C III, C II A1036.34, C II A1037.02, N III A991.58 


B.2 




360°-361° 


O 


VI, 


O V], Lya, II A1037.02, N III A991.58 


B.3 




350°-352° 


O 


VI, 


O V], Lya, C III, C II A1036.34, C II A1037.02, N III A991.58 


C.l 


15:07:18 


357°-358° 


o 


VI, 


C III 


C.2 




350°-352° 


o 


VI, 


O V], Lya, C III, N III A 989.90, N III A991.58 


C.3 




347° 


o 


VI, 


O V], [O V], Lya, C III, C 11 A1037.02, N 111 A991.58 


C.4 




346° 


o 


VI, 


O V], [O V], Lya, C III, C II A1037.02, N III A991.58 


C.5 




345° 


o 


VI, 


O V], [O V], Lya, C III, C II A1037.02, N III A991.58 


C.6 




344° 


o 


VI, 


O V], [O V], Lya, C III, C II A1037.02, N III A991.58 


D.l 


15:12:48 


347° 


o 


VI, 


O V], Lya, C 111, N III A991.58 


D.2 




346° 


o 


VI, 


O V], [O V], Lya, C III, C 11 A1036.34, C 11 A1037.02, N 111 A991.58 


D.3 




345° 


o 


VI, 


O V], [O V], Lya, C III, C II A1036.34, C II A1037.02, N III A991.58 


D.4 




344° 


o 


VI, 


O V], [O V], Lya, C III, N III A991.58 



Table 2. Line intensities of the O VI, O V, and electron density 



Blob No. 






O VI 








V] 


ne(xl0'5) 


Ne(xlOiS) 




Il032 


Vlos 


Vl/e 




Depth 


I1032/I1037 


I1218 


I1213/I1218 


V 


LASCO 


LASCO 


A 


34.0* 

39.4* x^=0.9 


-286 
-459 


58 
97 


33 
91 


0.858 


2.75±0.14 








2.04b 


13.6 


B.l 


261.6* 

132.0* x^=3.7 


-286 
-459 


81 
112 


63 
121 


1.217 


2.43±0.15 








1.31'' 


12.4 


B.2 


307.7 x^=0.9 


-487 


58 


33 


0.367 


2.78±0.19 


70.8±7.0 






5.14'' 


14.6 


B.3 


129.1 x^=2.8 


-260 


90 


78 


0.566 


2.05±0.15 


37.5±5.1 






2.38'' 


14.7 


C.l 


81.3 x^=0.2 


-457 


56 


30 


0.394 


3.03±0.32 








3.21'' 


9.68 


C.2 


118.8 x^=3.1 


-346 


57 


31 


0.401 


1.86it0.13 


35.6it4.9 






5.59'' 


17.1 


C.3 


574.7 x^=4.8 


-263 


71 


49 


0.501 


2.53±0.09 


62.5±6.5 


0.57±0.11 


1 28+°*^ t 


5.88'' 


22.4 


C.4 


621.4 x^=14.8 


-234 


87 


73 


0.613 


2.03±0.07 


39.5±5.2 


0.57±0.21 


1 25+l-''8t 
-^•^■^-0.53 


5.23'' 


24.3 


C.5 


71.1* 


-66 


38 


14 


0.652 


1.98±0.07 


49.5±5.8 


0.47±0.12 




5.73 


28.3 




439.2* x^=3.6 


-239 


54 


28 














C.6 


100.8* 


-95 


71 


49 


0.760 


2.05±0.10 


95.2±8.1 


0.10±0.04 


„ 9+10. 9tl) 


3.67 


21.1 




184.0* x^=2-0 


-239 


36 


13 














D.l 


176.3 x^=7.5 


-350 


68 


45 


0.537 


2.34±0.14 


54.7±6.1 






3.70'' 


15.0 


D.2 


167.3 x^=10.8 


-321 


48 


22 


0.378 


1.83±0.14 


43.3±5.4 


0.63±0.17 


1 06+°" t 


6.39'' 


18.2 


D.3 


142.0 x^=4.5 


-297 


30 


8.7 


0.234 


2.79±0.22 


86.7±7.7 


0.30±0.06 


3-08_o.72 


1.23 


19.9 


D.4 


82.3 x^ =3.8 


-297 


45 


20 


0.350 


3.43±0.42 


52.5±6.0 


0.37±0.10 




7.35'' 


19.4 



Note. — 

1 : Intensity (10* photons/cm^ sec sr) 
Vios, Vi/e : km / sec 
Tfc : lO^K 

Depth : Line of sight depth in Rq unit evaluated by (FWHM of the O VI) x (travel time from the source region to 

UVCS slit) 

* ; Two Gaussian fit 

t : The lower limit of electron density evaluated by [O V]/0 V] 
'' : The adapted electron number density for heating model 
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Tablc 3. Line intensities of the Lya, C III, C II, and N III 



Blob No. 


Lya 






C III 




C II (Int.) 


N III (Int.) 




^Lyot 




Vl/e 




Icj/J 




Vl/e 


1036.34 A 1037.02 A 


989.80 A 


991.58 A 


A 


4351.8 x^=222.4 


-380 


90 


4.9 


142.6" ±5.5 






2.6±0.7 3.2±0.8 




48.9= ±3.2 


B.l 


1776.9 x^=72.7 


-425 


64 


2.5 


334.3°±8.5 






7.7±1.3 20.5±2.1 




42.0±1.4 


B.2 


786.7±23.2 








a 






9.0±1.4 




48.4±3.2 


B.3 


5649.5 x^=51.8 


-325 


80 


3.9 


397.6±9.3 






5.8±1.1 20.5±2.1 




23.8±2.3 


C.l 










14.5±1.8 












C.2 


2169.3 x^=4.1 


-393 


70 


3.0 


285.2''±7.8 








6.3±1.2 


20.5±2.1 


C.3 


1226.5* 


-344 


56 


1.9 


220.5* 


-132 


30 


12.6±1.6 




26.0±2.4 




487.3* x^=34.1 


-231 


35 


0.74 


195.6* x^=4.0 


-193 


b 








C.4 


1254.3* 


-321 


67 


2.7 


40.4* 


-41 


b 


9.3±1.4 




18.1±2.0 




891.9* x^=22.0 


-208 


49 


1.5 


239.1* x^=7.4 


-193 


17 








C.5 


337.5* 


-317 


41 


1.0 


72.6* 


-46 


18 


12.2±1.6 




29.3±2.5 




1291.2* x^=31.9 


-136 


45 


1.2 


226.1* x^=3.7 


-198 


b 








C.6 


425.9* 


-317 


48 


1.4 


75.4* 


-76 


31 


5.0±1.0 




21.4±2.1 




565.0* x^=30.9 


-114 


35 


0.74 


190.4* x^=ll-4 


-198 


b 








D.l 


1442.6 x^=9-6 


-366 


39 


0.92 


285.8±7.9 










21.8±2.2 


D.2 


1441.0 x^=ll-6 


-344 


37 


0.83 


178.7±6.2 






3.8±0.9 6.6±1.2 




37.1±2.6 


D.3 


2762.1 x^=34.6 


-340 


41 


1.0 


312.0±8.2 






4.3±1.0 6.3±1.2 




35.3±2.8 


D.4 


1886.6 x^=19.5 


-317 


39 


0.92 


142.2±5.5 










16.2±1.9 



Note. — 

I : Intensity (10* photons/cm^ sec sr) 
Vios, Vl/e ■ km / sec 
Tfc : lO^K 

* : Two Gaussian fit 

" Observations at the edge of the observed wavelength range. Lower hmit of the C III emission. 
FWHM is within the instrumental width. Possibilities of " 
possible to be part of the Lya continuum. 
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Table 4. Energy budgets for the three heating models. 



Blob No. O VI (HI) Lya, C III, CII, N III (H2) 

QocQahh Qocn Qo^Qkr Qc^Qahh Qocn 





H.E. 


T.E. 


H.E. 


T.E. 




H.E. 


T.E. 


h 


H.E. 


T.E. 


H.E. 


T.E. 


A 


4.4-444 


0.1-9.7 


5.1-51 


0.1-0. 


,9 


87-106 


0.5 


81-84 


11 


0.3 


8.0 '' 


0.1 '' 


B.l 


129-420 


3.0-9.4 


57-90 


1.0-1. 


,6 


71-72 


0.3 


77 


13-21" 


0.3-0.5 " 


14-27" 


0.3-0.4" 


B.2 


58-65 


1.3 


57 


1.0 




71-77 


0.3 


77-78 


5.5-5.8 


0.1 


9.0 


0.2 


B.3 


15-89 


0.3-2.0 


9.0-57 


0.2-1. 


,0 


71-77 


0.3 


77-78 


5.2-5.5 


O.l*"* 


9.0 = 


0.2 = 


C.l 


18-225 


0.4-5.0 


25-100 


0.4-1, 


,7 


57-102 


0.2-0.4 


71-82 


2.5-45 


0.07-1.0 


4.0-25 


0.07-0.4 


C.2 


4.4-44 


0.1-1.0 


6.3-40 


0.1-0, 


,7 


57-70 


0.2-0.3 


71-75 


3.7-4.1 


0.1 


6.3 


0.1 


C.3 


60-68 


1.4 


63 


1.1 




62-70 


0.2-0.3 


73-75 


3.6-4.3 


0.1 


6.3 


0.1 


C.4 


74 


1.6 


63 


1.1 




70<= 


0.3^= 


75= 


4.2-4.7'' 


0.1'' 


6.3'' 


0.1'' 


C.5 


138-144 


3.4 


63 


1.1 




57-82 


0.2-0.3 


71-78 


5.4-5.7 


0.1 


6.3 


0.1 


C.6 


18-45 


0..3-0.8 


16-40 


0.3-0, 


,7 


57-82 


0.2-0.3 


71-78 


4.0-4.5'' 


0.09'' 


4.0'' 


0.07'' 


D.l 


51-56 


1.1 


17-44 


0.3-0, 


,7 


58-71 


0.2 


70-74 


4.9-5.3 


0.1 


6.9 


0.1 


D.2 


9.0-57 


0.2-1.0 


6.9-44 


0.1-0, 


,7 


58-105 


0.2-0.3 


70-81 


4.7-5.2'''= 


0.1'''= 


6.9''= 


0.1''= 


D.3 


6.7-65 


0.1-1.3 


6.9-44 


0.1-0, 


,7 


58-59 


0.2 


70 


6.2-6.3'' 


0.1'' 


6.9'' 


0.1'' 


D.4 


5.1-32 


0.09-0.5 


4.4-27 0.07-0.5 58-218 = 


0.2-0.7 = 


70-90= 


4.0-4.6= 


0.09= 


4.4 


0.07 



Note. — 

T.E. and H.E. ; Thermal Energy and Heating Energy (1.x 10^'* erg/g) assuming 10% Helium 
h : the percentage of the lost magnetic energy appearing as heat 
" : Does not satisfy the C II observations 

: Does not satisfy the N III observations 
= : Does not satisfy the O VI / O V constraint 

(and =) : Does not satisfy the C II < 3.* C II (and C III) criterion 
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C2 16:30:06UT 
Eir 1+:24:1 1UT 

UVCS 1+:17:56UT-19:16:05UT 




Fig. 1.— Composite image of EIT, UVCS and LASCO C2. The UVCS slit image shows 
the intensity distribution of the O VI hne summed over exposures from 14:17:56UT to 
19:16:05UT at 2.4 and P.A.=349°. The arrow indicates the active region where the 
bright flare was observed. 
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Fig. 2. — TRACE 1600 A at 14:27:07 UT. The arrow indicates the filament eruption. The 
lower part of the active region was seen as dark due to a saturation. 
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Fig. 3. — Observation time of each instrument. The arrows represent the first detection of 
this event at each instrument. In the TRACE 1550 A observation, the bright emission was 
already started at 14:04 UT. 
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+900 A 
91^ 88.9 86.1 



39.3 36.8 34.2 31.6 29.1 
+1000 A 






12.5 15.1 17.7 9.6 7.9 6.2 

+1200 A +970 A 
+900 A 

91.8 88.9 86.1 




59.3 36.8 34.2 31.6 29.1 
+1000 A 




12.5 15.1 17.7 9.6 7.9 6.2 

+1200 A +970 A 
+900 A 

91.8 88.9 86.1 




39.3 36.8 34.2 31.6 29.1 
+1000 A 




12.5 15.1 17.7 9.6 7.9 6.2 

+1200 A +970 A 

+900 A 
91.8 88.9 86.1 



Ji 



0.(1-4 




39.3 36.8 34.2 31.6 29.1 
+1000 A 



12.5 15.1 17.7 9.6 7.9 6.2 
+1200 A +970 A 



Fig. 4.— UVCS observations of the December 13 2001 CME. The left panel shows the O VI 
1032 A and 1037 A. The middle panel shows the Lya, [O V], O V], N III lines in both primary 
and redundant wavelengths. The primary (991.8 A - 986.1 A ) and redundant wavelengths 
(1212.5 A - 1217.7 A ) are represented in the upper and bottom axes in the middle panel, 
respectively. The right panel shows the C III line. . 
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Fig. 5. — Three-dimensional distribution of the O VI (red), Lya (green), and C III (blue) 
emission. Each grid presents 0.1 solar radius and the Sun is placed at the coordinate, (x:200, 
y:200, z:0). This figure is also available as an mpeg animation. The movie consists of two 
rotations, the first is the expansion of CME from 14:56 UT to 19:16 UT and the second 
rotation is at 19:16 UT. 
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Fig. 6. — Expansion heights of the CME plasma. The stars and crosses for C3 represent the 
position of the upper edge of left loop (where blob 'A' is placed) and right loop(where blob 
'C is placed) in the LASCO image in Figure 9. For the UVCS observation, the diamond 
represents blob 'A', the triangle represents blob 'B', the square represents blob 'C, and the 
mark, 'x', represents for blob 'D'. The solid hne is the distance applied for blob A, the dotted 
line is the distance apphed for blob B, and the dashed hne is the distance for blobs C and 
D. 
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Fig. 7.— Lya, [O V], O V], and N III line intensities of the blobs, C.5, C.6, and D.3. 
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Fig. 8. — The observed line ratio of the [O V]/0 V] and the electron number density evaluated 
from CHIANTI 5.2. 




Fig. 9 



, — The reconstructed positions of the 14 blobs are presented in the LASCO C3 image. 
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Fig. 10. — Temperature evolution of the protons and electrons of the model Qahh for 
matching the O VI and O V] line intensities of the blob C.5. 
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Fig. 11. — The line intensities when different heating rates are apphed to the heating model, 
QocQahh- The rainbow colors indicate different densities. The horizontal lines are the 
observed values. In the first panel, the solid and dotted lines are the observed O VI and O 
V values, respectively. In the upper third panel, the sohd and dotted lines are for the C III 
and C II lines, respectively. The green vertical line is the temperature that match the Lya, 
C III, C II, and N III emissions and the red vertical line is the temperature that match the 
O VI/ O V and O VI emissions. 
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Fig. 12. — The evolution of the density, temperature, and ionization states for one of the 
acceptable model for blob C.5. 
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Fig. 13. — Energy budgets for HI, Q,(xQ,ahh- 
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Fig. 14. — Heating Energy, o and + are for HI and H2, respectively. 
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Fig. 15. — Energy balance of the Kumar and Rust's heating model. 
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Fig. 16. — Time scale of radiative cooling, adiabatic cooling and thermal conduction. 



